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Abstract 

Exchange bias training effect in ferromagnetic/antiferromagnetic bilayers is investigated. In 
some systems the evolution of the exchange bias field He with the number of cycle n cannot be 
fitted by the empirical 1/y/n function. A unified expression is derived from a discretized Landau- 
Khalatnikov equation in the framework of the thermodynamics model which is proposed by Ch. 
Binek. This generalized model describes well training effect independent of the magnetization 
reversal mechanism in the ferromagnetic layers. 

PACS numbers: 75.30.Et; 75.30. Gw; 75.60.Jk 
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Exchange bias (EB) effect has been investigated extensively since its first discov- 
ery in CoO/Co nanoparticles jl-5]. After the EB is established in antiferro mag- 
netic (AFM) /ferromagnetic (FM) bilayers by (magnetic) field cooling from high temperature 
to below the Neel temperature of the AFM layer or by magnetic field during film deposi- 
tion, a shift of the hysteresis loop along the magnetic cooling field axis with an exchange 
field (He) is observed, companied by an enhancement of the coercivity (He). When cycling 
the bilayer through consecutive hysteresis loops, He and He often decrease monotonically 
with increasing cycling number n [Gjj. This so-called training effect has attracted extensive 
investigation because of its importance in both basic research and applications of spintronic 



devices 



2l|. 



During the training effect, large changes in He and Hq are often observed between n — 1 
and n = 2, which is thought to be strongly related to transition of AFM spins between easy 
axis [ll|. For n > 2, He decreases gradually with increasing n, which can be described by 



the empirical formula 
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where k is a parameter depending on the physics properties of the AFM layer. The empirical 
law was observed first for Co/CoO, NiFe/NiFeMn, and NiFe/C^Oa systems 6j. The EB 
training effect was considered in the framework of nonequilibrium thermodynamics by Ch. 
Binek, where the free energy of the system controls the relaxation process of the AFM 
spins towards the equilibrium state. With a discretized Landau-Khalatnikov equation, the 
following equation is obtained flo| . 



H E (n) - H E (n + 1) = 7(^(71) - H E (oo)) s 



(2) 



where 7 is the characteristic decay rate of the training behavior. Since Eq. [T] does not hold 
for some EB systems, such as NiFe/IrMn bilayers and Pt/Co/Pt/IrMn multilayers [tgI. \v\, 
several approaches have been proposed. For example, higher orders in the free energy were 
considered 15]. Alternatively, based on a thermal fluctuation model, a power law function 
on the cycle number was used to fit the training effect of He Moreover, the training 



effect was suggested to be strongly related to the evolution of the A 
training effect of He can be expressed by an exponential function 



? M spin disorder and the 
14, [La]. In this work, we 



will derive a unified equation to describe the variation of H E with n based on the theoretical 



model proposed by Binek, which can describe the EB training effect in FM/AFM bilayers 
with various FM magnetization reversal mechanisms. 

Two sets of samples were prepared by DC magnetron sputtering at ambient tempera- 
ture: glass/Cu(20 nm)/Ni 80 Fe 2 o(=NiFe)(6 nm)/Fe 5 oMn 5 o(=FeMn)(ll nm)/Au(10 nm) and 
glass/Cu(20 nm)/Fe a; Cr 1 _ ;c /Ir 2 5Mn75(=IrMn)(3 nm)/Cu(10 nm) with varying x. The base 
pressure of the system was 1.8 x 10~ 5 Pa and the working Argon pressure was 0.5 Pa during 
deposition. The deposition rates of the constituent layers were in the order of 0.1 nm/second. 
A magnetic field of 170 Oe was applied in the film plane during the film deposition to set 
the longitudinal exchange bias. In-plane magnetization loops were measured at room tem- 
perature by a Lakeshore vector vibrating sample magnetometer (VVSM). During training 
effect measurements, consecutive hysteresis loops were measured after deposition without 
break. All measurements were performed at room temperature. 

Figure [1] shows typical hysteresis loops at room temperature along the easy axis of 
NiFe/FeMn, Feo.51Cro.49/IrMn, and Feo.36Cr .64/IrMn bilayers. For the three samples, He 
at n = 1 equals 54.7, 14.5, and 24.5 (Oe) and is reduced by 28.9, 14.2, and 14.65 (Oe) after 
measurements of 160, 80, and 80 consecutive hysteresis loops, respectively. Large training 
effect can be found in all samples. Moreover, for NiFe/FeMn and Feo.51Cro.49/IrMn bilayers, 
the coercive field of the descending branch changes strongly while that of the ascending 
branch is hardly modified, resulting in a decrease of the coercivity. For Feo.36Cr .64/hMn 
bilayers, however, both branches shift towards the positive magnetic field direction and the 
coercivity almost does not change, i.e., from 28 to 25.5 (Oe) after 80 cycles (Fig. HJc)), as 
already observed for the perpendicularly magnetized Pt/Co/Pt/IrMn bilayers 17] . 

As shown in Fig. |2j the measured He(ti) in NiFe/FeMn, Feo.51Cro.49/IrMn, and 
Feo.36Cr .64/IrMn bilayers are fitted by using Eq. |2] One can easily find that the exper- 
imental results can be fitted by Eq. |2] Therefore, the EB training effect in these three 
samples can be described a discretized Landau-Khalatnikov equation. In principle, 7 and 
He(oo) can be fitted by using Eq. El Since He{ti) — He{u + 1) does not change much with 
n except for n = 1, however, the fitted values of these parameters are not rigorous. 

In order to examine whether the measured HeIii) can be described by the empirical law 
in Eq. [TJ H E for NiFe/FeMn and Feo.51Cro.49/IrMn bilayers is plotted as a function of l/\/n 
in Fig. |3] For different systems, 80 to 160 hysteresis loop cycles were measured to get an 
overall feature of the training effect. Only few tens of magnetization loops are usually mea- 
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sured 



Bflfl 



lOl . ll9M2l|. which is not adequate to inquire into the analytical expression of 
the training effect. In comparison, the data of epitaxially grown Fe/CoO bilayers are also 
given [if]]. Here, all samples exhibit the conventional training effect in which He decreases 
sharply between n = 1 and n = 2 and then changes slowly for large n JtJ. It is found that the 
measured data cannot be described very well by the empirical linear function of 1/ \fn even 
for n > 2. Based on these observations, we re-derived the dependence of He on the cycle 
number. Starting with Eq. [2] [lo|, the differential of He between neighboring consecutive 
hysteresis loops He{u) — He{ti + 1) is proportional to (He{u) — He(oo)) 3 when He changes 
slowly with n. Accordingly, one has the following differential equation, 

= - l{ H E {n) - H E (oo)f (3) 

where dn = 1 can then be obtained. Since He often changes slowly for n > 2, the analytical 
solution of the n-dependent He{ti) can be achieved as follows. 

H E {n) = He(oo) + \ (4) 



where n = ^h^^h^)) > strongly reflecting the fraction of the H E change between 
n = 2 and 3 in the entire training effect for n > 2. Similar parameters have been defined 
before j^]. If no = 2.0 in above equation, Eq. H] turns to the conventional empirical law, 
i.e., H E {n) — H E (oo) = i 2in - H E (oo), 7, and n can be fitted with the measured H E {n). 
Figures [3]^d)- [3^f) show that the fitting results and experiment data are in good agreement 
for NiFe/FeMn, Feo.51Cro.49/IrMn, and Fe/CoO bilayers. The value of no is fitted to be 6.6, 
5.8, and 0.97, and 7 to be 3.1 x I0~ 4 , 5.5 x 10~ 4 , and 7.0 x 10~ 4 (Oe -2 ), respectively. 

To further confirm its universal validity, Eq. H] is also used to analyze the training 
effect for Fen^fiCro^/IrMn bilayers and perpendicularly exchange-biased Pt/Co/Pt/IrMn 
multilayers 18|. Here, both samples exhibit the anomalous training effect 7j. Note that 
HeIu) of neither sample can be fitted with the 1/y/n function, as shown in Figs. HJ^a) 
& ID^b). Actually, for FeCr/IrMn bilayers the measured data at any n regime cannot be 
fitted by 1/y/n. However, the measured data can be described by Eq. H]very well as shown 
in Figs. ET|c) &ET|d). The n is found to be 10.3 and 6.5 and 7 is 2.3 x 10" 4 and 1.3 x 10~ 6 
(Oe -2 ) for Fe/CoO and Pt/Co/Pt/IrMn multilayers, respectively, it is therefore evident 
that the model based on the AFM spin relaxation process holds for training effects in 
various FM/AFM bilayers with either the conventional training effect or the anomalous 



training effect. 

When the magnetization reversal is dominated by domain wall depinning, such as in 
Feo.36Cr .64/IrMn and Pt/Co/Pt/IrMn systems, He shows the so-called anomalous training 
effect in which He(1) — He (2) is small, i.e. the first-jump vanishes fll |. Meanwhile, when 
the magnetization reversal is dominated by magnetization coherent rotation, such as in 
NiFe/FeMn, Feo.51Cro.49/IrMn, and Fe/CoO bilayers, He exhibits conventional training 



effect in which the first-jump appears 



ff 



16-18 



|23| . Therefore, the magnetization reversal 



mechanism plays an important role in the training effect. Moreover, when the evolution of 
H E {n) deviates seriously from the l/\/n linear dependence, such as for Pt/Co/Pt/IrMn, 
FeCr/IrMn, and NiFe/FeMn systems, a large no is acquired. Otherwise, a small no is 
obtained, such as for the Fe/CoO bilayer. Furthermore, according to n = 2 (H^2)-^i(3)) 
and Eq. [2], one has H E {2) — H E (oo) = , that is to say, the n 7 product determines the 
training effect for n > 2, i.e., H E {2) — H E (oo). Finally, with the values of the parameters 7 
and no, the n dependence of H E {n) can be identified. 

In summary, for NiFe/FeMn, Feo.36Cro.64/IrMn, and Feo.51Cro.49/IrM bilayers the magne- 
tization reversal process is dominated by domain rotation, the domain wall depinning, and a 
combination of both modes, respectively. Depending on the reversal process, the evolution 
of H E {n) may be close to or seriously deviate from the empirical l/^/n function. Based 
on the thermodynamics model proposed by Ch. Binek [10], a unified analytical function 



of l/\/n + uq — 2 is derived for n > 2 and is verified to be applicable to the thermally 
activated EB training for various FM/AFM bilayers. 
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FIGURE CAPTIONS 



Figure 1 Typical hysteresis loops of NiFe/FeMn (a), Feo.51Cro.49/IrMn (b), and Feo.36Cro.64/IrMn 
(c) bilayers. The inset numbers refer to the cycle numbers of hysteresis loops. 
Figure 2 The curve of the measured He{ti) — He(ti + 1) versus He(ti) for NiFe/FeMn (a), 
Feo.51Cro.49/IrMn (b), and Feo.36Cro.64/IrMn (c) bilayers. The solid lines refer to the fitted 
results based on Eq. [2j 

Figure 3 The n dependence of the measured He (o pen symbols) for NiFe/FeMn (a,d), 



171 ] (c,f) bilayers. The measured data were fitted 



Feo.51Cro.49/IrMn (b,e), and epitaxial Fe/CoO 
using a linear function of l/y/n (the left column) and l/y/n + riQ — 2 (the right column), 
respectively. The solid lines refer to the fitted results based on Eq. [H 



Figure 4 The n dependence of the measured He for Feo.36Cro.64/IrMn bilayers (a,c) |16 : ] and 
Pt/Co/Pt/IrMn multilayer (b,d). The measured data were fitted using a linear function of 1/y/n 



(the left column) and l/\/n + no — 2 (the right column), respectively. The solid lines refer to the 
fitted results based on Eq. fH 
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